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61,   Introduction 

In  many  cases  of  practical  importance  there  arises 
the  problem  of  determining  the  asymptotic  behavior  for 
large  values  of  n  of  a  complex  Integral  of  the  form 

(1.1)  F(n,a)  =  Je-^^^''^^  i(z)dz 

C 

where  n  is  a  positive  real  nximber,  f  and  ^  are  meromorphic 
functions  of  the  complex  variable  z,  f  is  continuous  in  the 
complex  parameter  a,  $^(z)  is  of  at  most  exponential  growth 
at  infinity  and  C  is  a  contour  in  the  z-plane  which  extends 
to  infinity  in  both  directions  and  on  which  Re  f  >  0 
sufficiently  far  out  (to  ensure  convergence  of  the  integral), 
Some  of  the  assumptions  about  n  and  C  are  not  essentialf 

If  n  is  complex,  fut  has  a  fixed  argument  9,  we  only  have 

I  I  is  18 

tp  write  n  =  |n|e*  and  absorb  the  e   into  the  function  f. 

If  Q   has  one  or  both  end  points  at  finite  points,  then  it 

Vill  be  seen  that  the  following  discussion  need  undergo 

only  trivial  alterations, 

A  method  of  attacking  the  problem  has  long  been  known 

and  often  utilized.   We  shall  discuss  this  method  in  greater 

detail  in  the  next  section.   In  its  simplest  form,  the  method, 

often  called  "the  saddle  point  method",  consists  in  first 

making  the  change  of  variable 

(1.2)  K  =   f (z,a) 

in  (1.1)  and  then  showing  that  the  asymptotic  behavior 
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can  be  obtained  by  finding  the  branch  points  z.  of  the 
mapping,  i.e.  the  points  where  f  (zv^,a)  =  0,  taking  their 
images,  ^.  ,  in  the  ^-plane,  expanding  z  and  therewith 
|(z)4p  into  series  in  fractional  powers  of  (^-^^),  and 
Integrating  termwise.   The  results,  as  will  be  shown  in 
the  next  section,  are  series  of  the  form  e     times  a 
power  series  in  powers  of  n    where  v  is  the  order  of  the 
first  non-vanishing  derivative  of  f  at  z.  .   Due  to  the 
presence  of  the  exponential  factor,  it  is  clear  that  the 
branch  point  ^,  with  the  smallest  real  part  will  contribute 
the  dominant  series  to  the  asymptotic  expansion.   (It  will 
be  seen  later,  however,  that  only  branch  points  ^^  lying  in 
the  ^-plane  to  the  right  of  the  image  curve  I  of  C  con- 
tribute anything  to  the  expansion. ) 

There  is  one  serious  difficulty  in  the  method, 
however.   As  alluded  to  above,  the  nature  of  the  series 
for  z  in  the  neighborhood  of  i^^   depends  upon  the  order  of 
the  first  non-vanishing  derivative  of  f  at  z^.      If,  as 
indicated  by  the  notation  in  (1.1)  and  (1.2),  f  happens 
to  depend  upon  another  parameter,  a,  then  the  order  of  this 
"first  non-vanishing  derivative"  at  a  given  branch  point 
may  be  different  for  different  values  of  a.   From  the  point 
of  view  of  practical  applications  the  most  Interesting 
case  Is  that  in  which  the  order,  v,  is  the  same  for  all 
values  of  a  in  a  neighborhood  of  a  particular  value,  say 
a  =  0,  but  becomes  greater  than  v  when  a  =  0;  i.e. 
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f^v)(z,a)  ^  0  for  all  a  ?^  0,  but  f(^^z,0)  =  0.  I  ^t^^  ^^^  ' 

For  the  sake  of  simplicity  and  definiteness  we  shall 
treat,  in  the  present  paper,  the  case  in  which  the  branch 
point  is  of  first  order  (v=2)  when  a  ^  0  and  of  second 
order  (v=3)  when  a  =  0,   More  complicated  cases  seldom 
arise  in  practice  and,  in  any  event,  it  seems  unlikely 
that  branch  points  of  higher  order  could  cause  any  real 
theoretical  difficulty,  for  the  idea  of  the  method  to  be 
presented  seems  capable  of  straightforward  generalization. 

Thus,  let  f(z,a)  be  a  function  analytic  in  z  in  a 
neighborhood  of  the  origin,  z  =  0.   Let  z_|^(a)  be  the  two 
(distinct)  zeros  of  f'(z,a)  closest  to  the  origin.   (Here 
and  henceforth  primes  denote  differentiation  with  respect 
to  the  variable  other  than  a  and  we  assume  that  for  a  7^  0 
there  are  exactly  two  such  zeros);  i.e.  let 

(1.3)       f'[z^(a),a]  =  0 

where  2^(a)  ^   z_(a). 

Assume,  however,  that 

(l.ij.)      z^(0)  =  z_(0)  =  Zq  (say), 

so  that  for  a.  ^  0 , 

(1.5)  f"[z^(a),a]  ^   0  , 
but,  at  a  =  0, 

(1.6)  f"[z^(0),o]  =  0. 
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We  require,  though,  that 
(1.7)     f'"[z^(a),a]  ^  0      for  all  |a|i  <  A 

(Including  a  =  0). 

In  this  case,  apart  from  an  exponential  factor,  the 
saddle  point  method  yields  an  asymptotic  expansion  in 
powers  of  n  ^  for  a  7^  0,  but  one  in  powers  of  n   -^  if 
a  =  0.   The  coefficients  in  the  expansion  for  a  j^  0  are 
naturally  functions  of  a  and  one  might  well  ask  what  happens 
to  them  if  we  let  a  —s.  0.   The  answer  is  not  too  surprising 
in  view  of  the  fact  that  the  series  for  a  =  0  proceeds  in 
different  powers  of  n  from  that  for  a  7^  0;  the  answer  is 
that  they  blow  up.   In  fact,  the  three  different  expansions, 
one  for  a  <  0,  one  for  a  =  0,  and  one  for  a  >  0  in  general 
have  to  be  computed  separately. 

Now  it  is  well  known  that  one  of  the  principal  uses 
of  asymptotic  series  is  for  the  purpose  of  numerical  compu- 
tation and  it  might  be  expected  that  the  existence  of  three 
different  asymptotic  series  in  n  for  three  different 
ranges  of  the  parameter  a  would  lead  to  poor  numerical 
results  for  values  of  a  very  close  to  a  =  0,   This  is  indeed 
the  case.  As  a  result,  numerous  efforts  have  been  made  to 
procure  asymptotic  representations  valid  in  the  so-called 
"transition  region"  when  a  is  close  to  zero.   Most  of  these 
efforts,  however,  have  been  directed  toward  solving 
specific  problems  rather  than  the  general  one  and,  as  a 
result,  their  derivation  generally  requires  the  use  of 
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highly  specialized  techniques,  techniques  which  can  rarely 
be  extended  to  many  examples  outside  the  particular  case 
for  which  they  were  developed.   Moreover,  they  seldom  if 
ever  yield  full  asymptotic  expansions,  but  rather  only 
one  term  or  so,  and  their  asymptotic  character  is  often 
difficult  to  prove.   It  would  therefore  be  highly  desirable  j 
to  have  a  uniform  (in  a)  asymptotic  expansion  valid  for  all 
a  in  a  neighborhood  of  a  =  0. 

Only  very  recently,  F.  w.  j.  oiver  [9]'"  has  succeeded 
in  obtaining  such  uniform  expansions  for  those  integrals 
of  the  type  (1.1)  which  are  solutions  of  a  certain  class 
of  second  order  linear  ordinary  differential  equations. 
Actually,  he  is  interested  in  expansions  for  solutions 
for  such  equations  and  not  for  integrals  of  the  type  (1.1) 
so  that  in  order  to  apply  his  results  to  an  integral  like 
(1.1)  it  would  first  be  necessary  to  find  the  differential 
equation  (if  it  exists  -  it  may  not)  satisfied  by  (1.1). 
In  some  cases,  such  as  Bessel  functions,  this  can  indeed  be 
done.   In  others,  it  cannot.   It  would  therefore  still  be 
quite  desirable  to  obtain  a  uniform  expansion  for  an 
integral  such  as  (1,1). 

In  the  present  paper  a  method  for  obtaining  such  a 
uniform  expansion  will  be  presented.   Differing  considerably 
from  Olver's  method,  it  consists  in  replacing  the  usual 
change  of  variable,  (1,2),  by  a  slightly  more  complicated 

Numbers  in  square  brackets  refer  to  the  bibliography  at 
the  end  of  the  report. 
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one,  namely 

(1.8)  b  +  3a^w  -  w^  =  -f(z,a) 

where  w  is  the  new  variable  of  integration  and  where  (the 
reasons  for  the  following  choice  of  a  and  b  are  explained 
at  length  in  §3. ) 

(1.9)  a^  =  J[f(z.,a)  -  f(z^,a)] 

(1.10)  b  =  -  i[f(z^,a)  +  f(z_,a)] 

(cf.  (1.3)  -  (1.7)  for  the  definition  of  z^). 

Under  the  change  of  variable  (1.8)  the  integral 
(1.1)  goes  over  into 

(1.11)  F(n,a)  =  e^^  ;  e'^^^^^^"^^^  cKw,a)dw  , 

K 

where 

(1.12)  4(w,a)  =  $[z(w,a)]  ^   , 

z(w,a)  being  defined  implicitly  by  (1.8),  and  where  K  is 
the  image  of  C  under  (1.8). 

As  a  result  of  the  choice  of  a  and  b  the  mapping  (1,8) 
can  be  (and  in  8i|  is)  proved  to  be  one  to  one  in  a  neigh- 
borhood including  both  branch  points  z_^_   provided  a  is 
sufficiently  small.   This  means  that  z  can  be  expanded  into 
a  convergent  power  series  in  w  in  this  neighborhood. 
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Slightly  more  than  this  is  needed.   It  is  shown  that  4  can 
be  expanded  into  a  convergent  series  of  the  following  type: 


00 


(1).  2  ,2xk 


c/.^^a^.w^)^]  , 


(1.13)    (t(w,a)  =  J2     i<   'i^   '^   )      +  w  ., 

k=0    ^  ^ 

the  coefficients  of  which  are  not  too  hard  to  compute; 
iie.  a  prescription  for  computing  the  coefficients  can  be 
given.   In  practice  it  might  be  rather  difficult  to  compute 
more  than  the  first  few,  but  theoretically  it  can  be  done. 

The  path  K  of  the  integral  in  (1.11)  is  then  deformed 
so  as  to  pass  through  the  region  of  validity  of  the  series 
(1.13)  and  after  a  number  of  remainder  estimates  and  some 
further  path  pushing  it  is  shown  that  termwise  integration 
of  the  series  (1,13)  over  the  path  L,  =  I^  -  I,  shown  below 

w-plane 
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yields  an  asymptotic  series  of  the  form 

(l.l!|)     F(n,a).>..  e""^  XI  C,(n,a) 

k=0  ^ 

where  the  c)    (i  =  1,2)  are  the  coefficients  appearing 
in  the  series  (1,13)  and  where 

(1.15)  A^(n,a)  =   f  (a^-w^)""  ^niZ^^^-^h   ^„ 
and 

(1.16)  Bj^(n,a)  =  j   w(a2-w2)^  ^nOaVw^)  ^^  ^ 

Note  that  A,  and  B,  are  really  functions  of  n  and  a  since, 
as  remarked  previously,  a  is  a  function  of  a. 

If  we  define  A_,  =  B_,  =  0,  then  the  integrals  A, 
and  Bj^  satisfy  recursion  relations  of  the  form 

(1.17)  Aj^(n,a)  =  -^i^  Bj^.^^*^'^^ 
and 

(1.18)  Bj^(n,a)  =  -  I2|2li  A^_^(n,a)  ^  ^^^^^   A^_2(n,a) 

for  k  >  1. 
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By  using  (1.1?)  and  (1.18)  it  will  be  shown  later 

that  the  functions  A,  and  B,  and  therewith  C^  (cf.  (l.li|)) 

are  of  decreasing  order  of  magnitude  in  n.   For  example, 

k 
/  C,  (n,a)  =  0(n  ^)   CQ(n,a)      if  k  is  even 


(1.19)   ',  1,^1 


(   C^(n,a)  =  0(n   "^    )   CQ(n,a)    if  k  is  odd 

independently  of  a  (see  p.  ^2),      This,  plus  an  estimation  of 
the  remainder,  (which  is  also  carried  out  later),  establishes 
the  asymptotic  character  of  the  series  (l.lij.). 

The  recurrence  relations  (1.17)  and  (1.18)  also 
enable  one  to  express  all  the  A,  '  s  and  B,  's  in  terms  of  A_ 
and  B_,   Of  course,  all  this  would  be  of  very  little  use  if 
Aj-.  and  B^-.  were  not  themselves  expressible  in  terms  of  some 
well  known  functions,  functions  sufficiently  well  tabulated 
to  enable  the  formula  (1.11+)  to  be  of  use  for  numerical 
computations.   Fortunately  (and  in  fact  this  is  one  of  the 
reasons  for  the  choice  of  the  transformation  (1.8))  A_  and 
Bq  are  expressible  in  terms  of  well  known  functions.   If 

zt-  W 


,     (         2t-  i 
,i  iz)    =   — i—   i   e     ^ 


(1.20)      Ai^(z)  =  ^    J       e  ^   dt 

where  L_j^  is  the  same  path  as  shown  on  p.yi  is   ome  of  the 
Airy  integrals  (see  p.  B-67  of  Friedrichs'  notes  on  Special 
Topics  in  Analysis  [I4.]  or  Jeffreys  and  Jeffreys  [?],  p.  ij.76), 
then 
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(1.21)   A_(n,a)  = 


and 


27ti 


0' 


(3n) 


n)V3    ^ 


AlJ(3ny^  &^] 


(1.22)   B^(n,a)  =  ^""^^^   Al|[(3n)^^  a^] 


0' 


(3n) 


^/3 


(prime  denotes  derivative  with  respect  to  the  argument). 

For  the  purposes  of  nvunerical  computation  it  is 
somewhat  more  convenient  to  relate  A»  and  B^  to  the 
following  integral  which  is  also  a  solution  of  Airy's 
differential  equations 

,    ^   zt+  it^ 
(1.23)   h^(z)  =1^/6    ^   dt 


(k  here  is  a  numerical  constant  and  has  no  relation  to 
any  indices  appearing  in  preceding  formulae)  where  L,  is 


the  path  shown  below 


t -plane 
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and  7-  <  w  <  -^  (again  this  w  has  no  relation  to  our  w). 


The  integral  h, (z)  together  with  its  derivative,  h, (z) 
is  tabulated  for  complex  values  of  the  argument  z  in  the 
Harvard  Computational  Laboratory's  "Hankel  Functions  of 
Order  /3  "  [5].  A_  and  B^  are  expressed  in  terms  of 
h,  and  h,  by  the  following  formulas: 

(1.2i^)    AQ{n,a)  =  "^  ^  .^y   h^[e  "^  (3n)^  e^]  ; 

k(3n)  ^^ 

(1.25)    BQ(n,a)  =     "^  ^  h^  [e'  ^(3n)^  a^] . 

k(3n)  6 

Using  all  these  formulas,  theoretically  any  desired 
number  of  terms  in  the  series  (l.li;)  can  be  written  down, 
although  in  practice  seldom  more  than  the  first  term  is 
needed.   For  it,  only  Cq   (a)  is  needed,  and  (see  §6) 


(1.26)  cil)(a)=|j(z,)/7=f|=;  .|^(._,  pSEZ; 

VZZ+  VZZ-' 

If  we  let  a  —?»  0,  than  a  — ?»  0  as  does  f   and  c^  '(a) 
goes  over  continuously  into 

3 


U.27)   o<i>,o,M,.„,.^,--^^ 

Thus,  by  (l.li^),  (1.2)4)  and  (1.26), 
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(1,28)   F(n,a)^ 


n   e 


i-rt 


2k(3n) 


V^ 


$(-.^/f 


TT 


zz 


(z+»a) 


V    22   -"'     J 


in 
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which  is  an  asymptotic  formula  valid  for  large  n  and  s-nall  a. 
In  s8  formula  (1,28)  is  applied  to  the  Kankel  functions 
of  nearly  equal  order  and  argument  to  yield  the  formula 


(1.29)   H^^^nr).-..- 


-i 


-6(.  r'^-l-tan  Vr'^-l)  ^^ 


k(3n) 


V3,./ 


y^yT^ 


*   h. 


(3n)^''^  (>/r^-l-tan-Vr^-i)^/3 


.^3 


in  which  (as  will  be  explained  later)  it  is  necessary  to  take 
the  above  form  for  r  >  1  and  to  write  v  r   -1  =  ii/l-r 
and  choose  i    =-iifr<l.   i-ormula  (1.29)  is  valid 
for  large  values  of  n  and  for  r  close  to  1, 
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§2.   The  Saddle  Point  Method 

As  mentioned  in  the  introduction,  there  is  a  well 
known  standard  way  for  obtaining  an  asymptotic  expansion 
for  an  integral  of  the  type  (1.1),   This  Tiethod,  which  was 
originated  by  Debye  [3]  in  1909,  yields,  apart  from  an 
exponential  factor,  an  asymptotic  series  for  F  in  fractional 
powers  of  n.   As  originally  conceived  by  Debye  and  later 
elaborated  upon  by  others  the  method  possessed  labyrinthal 
arguments  about  the  nature  of  so-called  "paths  of  steepest 
descent",  paths  passing  through  the  cols  of  Re  f  and 
on  which  Re  f  decreases  most  rapidly  upon  leaving  the 
cols.   In  addition  to  the  difficulties  in  understanding 
the  ideas  of  the  method,  there  are  formidable  practical 
difficulties  in  computing  the  coefficients  in  the  expansion. 
To  date  no  one  has  succeeded  in  obtaining  any  essential 
simplification  of  the  computational  difficulties.   However 
a  considerable  simplification  in  the  presentation  of  the 
ideas  of  the  method  has  been  achieved  recently  by  B.  L, 
van  der  Waerden  [10]  who  describes  the  saddle  point  method 
(as  it  is  often  called)  as  follows? 

In  (1.1)  introduce  the  new  complex  variable  ^  =  f(z,a) 
as  in  (1.2).   Then  C  goes  over  into  some  curve, j  ',  in  the 
w-plane  which,  because  of  the  assumption  that  Re  f  >  0 
far  out  on  C,  goes  to  infinity  in  the  right  half  of  the 
^-planej 


ve.v 


1.1 


lu. 


^-plane 


The  integral  (1.1)  goes  over  into 

(2.1)    F(n,a)  =  '   e"""^  \{z{K,^\   ||  d? 


or 


(2.2)  F(n,a)  =  ;   e"^^  >|/(^,a)  d?  , 
v;here 

(2.3)  ^(^,a)  =  ({[z(^,a)]  ||  . 

If  the  relatively  mild  assumption  that 

(2.i+)     k(^,a)|  <  Me^'^'    (M  and  k  constants) 

for  Re  Z,      sufficiently  large  is  made,  then  the  path  of 
integration,;  ',  can  be  deformed  in  such  a  way  as  to  yield 
the  desired  asymptotic  series  for  F.   It  is  well  knovm  that 
the  values  of  complex  integrals  are  determined  by  the 
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singularities  of  their  integrands.   Clearly  the  only  part 
of  the  integrand  of  (2.2)  that  can  have  any  singularities 
is  the  function  ij/(z,a)  and  it  is  obvious  from  (2,3)  and 
the  assumption  that  f  and  ^  are  meromorphic  that  the  only 
possible  singularities  of  \|/  are  the  poles  of  ^  and  the 
branch  points  of  the  mapping  (1,2),   f^irtherraore,  only 
those  poles  and  branch  points  which  lie  to  the  right  of  •'"* 
will  be  of  importance  in  evaluating  the  Integral  (2,2). 
For  each  branch  point  to  the  right  of  '  we  introduce  in 
the  ^-plane  a  cut  parallel  to  the  real  axis.   We  assume 
for  simplicity  that  no  such  branch  cut  meets  the  curve  /  , 
for  otherwise  we  could  always  modify  the  cuts  in  order 
to  achieve  this. 

Once  this  has  been  done  it  is  then  clear  by  Cauchy's 
theorem  that  the  integral  {2,2)    is  equivalent  to  a  sum  of 
integrals  over  the  banks  of  the  branch  cuts  plus  a  sum  of 
residues  from  the  poles, 

2;-plane       J^^ 
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If  ^^   is  a  typical  pole,  say  for  simplicity  of 
order  one,  then  it  contributes  a  terra  of  the  form 

(2.5)    Fp{n,a)  =  e   P  fp^^p*^) 

(p  here  is  merely  a  subscript;  no  differentiation  is  iiiplied) 
to  the  asymptotic  expansion  of  F, 

On  the  other  hand,  a  branch  point,  ^,  ,  contributes 
a  terra  of  the  form 

'b 

where  L,  is  a  loop  surrounding  the  branch  cut  from  ^,  to 


(2.6)    F^(n,a)  =  j   e"^^  »|»(?,a)d?  , 


infinity.   Now  f(z,a)  is  by  hypothesis  analytic  in  z. 
Again  for  siraplicity,  assume  that  f  is  analytic  at  the 
pre-image,  z,  ,  of  the  branch  point  ^,  .   Then 

(2.7)  ?  =  f(z,a)  =  CQ(a)  +  c^(a)(z-z^)^  +  ... 
where 

(2.8)  c^(a)  ^  0    and  v  ^  2. 

Put^^^^ 

/■2 


(2.9)  U  iz,a)   =  j    §(r)dr 


^b 


This  idea,  which  avoids  the  necessity  of  inverting  the 
series  for  \|/  and  ^  separately,  was  suggested  by  Prof, 
A.  S.  Peters.     °^ 


^  > 


J    :"'     ' 


17. 
Then 
(2.10)   c.;^(z,a)  =  4[z(^,a)]  ||  =  »|/(^,a)  . 

Now  by  the  La grange -Burmann  formula  [6],  the  series  (2.7) 

can  be  inverted  and  we  can  expand  0(7,0)  in  a  series  in 
1 

V 


(^-^b>   •• 


li. 
00  — 


(2.11)   >./.!  (z,a)  =  ^^  r^(<^)ii:'K^)       » 


[1=0 


whence 


CO  — 

(2.12)    i|/(2,a)  =  l2     r,A^)iK-Ky,)^ 

ix  =  0  ^  ° 


• 


The  series  (2.12)  will  converge  absolutely  and  uniformly  in 

some  neighborhood,  |^-?v|  <  B,  of  ^,  •   In  such  a  neighborhood 

we  can  insert  it  into  (2.6)  and  integrate  term  by  term, 

obtaining 

CO         '^  - 

(2.13)    F;(n,a)  =  ^I  Y^(a)   '^   e-"^(^-?j^)^  d^  , 

where  L,  is  that  part  of  the  path  L  which  lies  in  the 
neighborhood  in  question. 

Next,  in  order  to  be  able  to  evaluate  explicitly  the 
integrals  occurring  in  (2.13)  we  extend  the  path  L,  to 
infinity  again  (i.e.  replace  it  by  Lv^).   Such  an  extension 

is  clearly  permissible  in  virtue  of  the  presence  of  the 

-n^ 

exponential  e  ^  in  the  integrand.   Of  course,  the  resulting 
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:a)    )     e" 


(2.11;)   S  (n,a)  =  ^  y  (a)    )  e'^^C^-?,  )^  d^ 

.=0  ^     .^ 


is  no  longer  equal  to  F.  .   It  can  be  shovm,  however,  and 
indeed  we  shall  show  it  later  on  in  a  slightly  different 
context,  that  the  series  S.  is  an  as3n7iptotic  expansion 
for  F,  as  n  — «» co  •   Moreover,  the  terms  of  S,  are  easy  to 
compute,  for  on  the  part  of  L,  which  runs  along  the  upper 
bank  of  the  branch  cut  we  can  use  the  parameter 

1     1 

(2.15)      (^-^b)"  =  r"    ,    0  g  r<  00, 

whereas  on  the  lower  bank  we  can  take 

1    2l2t  1 

(2.16)  iK-K^)"  =  e  "    r"   ,      0  g  r<  00. 

The  series  (2.1i4.)  then  becomes 

(2.17)  S,  (n,a)=2:  Y^(a)(l-e    ^     )     [       e  ^  f''   d  T 


-n?v   00 


0 

2i7i'^  CO  a 


0 


=  e        -   ^   u-e    -      ;   Y,,U;     »        e  "•   l'-    at 


r.y        rr.  2iTt|J.  CO  ^  -  1 
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n 


Thus, 


-n^,  oo  ^     r„(a)r(^) 


b 


(2.18)      S^(n,a)   =e        ''^(l-e^) 


ii 
It    is   clear   from   (2.18)    that   the  branch  point  ^, 


lJ-=«  V 

n 


'b 
with  smallest  real  part  (i.e.  closest  to  D  vjill  make  the 

major  contribution  to  the  as^^miptotic  expansion  of  F,  for 

the  other  branch  points  will  all  contribute  series  which 

are  multiplied  by  factors  exponentially  smaller  than  that 

of  the  dominant  series.   If  ^  is  the  branch  point  with 

smallest  real  part,  then 

-n^  00      ^2lii  Y  (a)r(H) 
(2.19)   P(n,a)^-^e   ^  J"  (1  -  e  ^  )  -^i—    ^ 


^  n^ 

is  an  asymptotic  series  for  F(n,a).   Of  course,  our  above 
discussion  bv  no  means  constitutes  a  rigorous  proof  of 
this  fact,   ''''e  have  merely  tried  to  indicate  the  salient 
features  of  the  situation  because  below  we  shall  have  to 
give  a  detailed  proof  of  the  asymptotic  character  of  a 
somewhat  different  series  and  the  method  of  proof  will  be 
almost  the  exact  analogue  of  the  one  given  by  van  der 
Ifaerden  for  the  case  above. 
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83.  A  Difficulty  in  the  Saddle  Point  Method 

Having  completed  our  sketch  of  van  der  Vaerden's 
exposition  of  the  saddle  point  method,  we  come  nov;  to  the 
problem  which  is  the  reason  for  the  present  paper.   It 
will  be  noted  that  although  vie   have  assumed  that  the  in- 
tegral P(n,a)  given  by  (1.1)  depends  upon  a  parameter  a 
as  well  as  upon  n,  we  have  up  to  the  present  time  made  no 
ostensible  allowance  for  this  fact.   There  has,  hov-jever, 
been  one  crucial  point  where  the  parameter  a  entered  in  a 
non-trivial  way,  namely  in  the  assertion  (2.8).   Since 
c  (a)  is  continuous  in  a  and  ?^  0,  it  is  clear  that  there 
will  be  ranges  of  a  wherein  c  (a)  ^  0.   There  may,  however, 
be  some  value  of  a  for  which  c  (a)  =  0.   V'e  assiome  that 
this  value  of  a  is  isolated  and  without  loss  of  generality 
may  take  the  point  to  be  the  origin,  a  =  0,  i.e.  we  suppose 
that 

(3.1)  c^(0)  =  0  , 
but 

(3.2)  Cy(a)  /  0        for   0  <  |a|  <  A     (A  >  0). 

Then  it  is  not  hard  to  see  that  the  coefficients, 
Y  (a),  of  the  series  (2.12)  blow  up  as  a  — >  0.   In  other 
words,  the  entire  procedure  breaks  down.   Actually,  things 
are  not  quite  that  bad,  since  for  any  fixed  value  of  a  an 
asymptotic  expansion  can  always  be  obtained  in  the  form 
(2.19),  whereas  for  a  =  0  a  similar  expansion  with  different 
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-  JL 
v' 

coefficients  and  proceeding  in  powers  of  n     where  v* 

is  the  smallest  integer  >  1  such  that  c  , (0)  ^  0  (such 
an  integer  surely  exists)  can  be  found.   Hov/ever,  the 
expansions  for  a  ^  0  (there  are  in  general  two  of  them, 
one  for  a  >  0  and  a  different  one  for  a  <  0)  do  not  go 
over  continuously  into  that  for  a  =  0  as  a  —^  0  in  them. 
In  fact,  as  mentioned  previously,  they  do  not  go  over 
continuously  into  anything;  they  blow  up. 

As  explained  in  the  introduction,  there  is  need 
for  a  uniform  expansion.   I ith  a  view  toward  arriving  at 
such  an  expansion,  let  us  analyse  a  little  more  closely 
the  situation  described  by  (3.1)  and  (3.2).   From  (2.7) 
we  see  that 


dz 


'f(2^,a) 


Z=Zi 


kl 


■;z 


Relations  (3.1)  and  (3»2)  therefore  mean  that 


(3.1|) 


^  f(z.  ,a) 


\ 


t 


=   0 


£  =   1,2,...  v-1,  for  all  a. 


while 


(3.5) 


^  0 


if  a  7^  0. 


For  a  =  0,  however. 
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(3.6) 5 =  0  . 

0  z 

In  other  words,  z^  is  a  branch  point  of  order  v  for 
a  7^  0  but  becomes  a  branch  point  of  order  >  v  for  a  =  0, 
and  it  is  precisely  this  phenomenon  that  causes  the  break- 
down of  the  series  expansion  (2.12)  and  the  resulting 
failure  of  the  saddle  point  method. 

As  mentioned  in  the  introduction,  we  shall  treat 
only  the  case  in  which  two  branch  points  of  second  order 
coalesce  to  become  a  branch  point  of  third  order,  i.e.  in 
which  assumptions  (1.3)  -  (1.7)  are  satisfied.   From  them 
it  follovjs  that  in  the  neighborhood  of  the  origin 

2         3 

(3.7)  f(z,a)  =  CQ(a)  +  c^(a)z  +  C2(a)z  +  c,(a)z-^  +  ... 

where  c^(a),  C2(a)  ;^  0  for  a  7^  0  but  c^(0)  =  0^(0)  =  0, 
while  c^(a)  7^  0  for  all  a  in  |a|  <  A  (including  a  =  0). 
Thus,  if  we  put  ^  =  f(z,a),  then  z  is  a  pov;er  series  in  ^ 
for  a  ^t  0,  but  a  power  series  in  J^  3  if  a  =  0.   I-Tiat  is 
clearly  needed  is  a  uniform  expansion  for  z  in  powers  of  a 
variable,  say  w,  an  expansion  which  retains  its  validity 
even  for  a  =  0.   Obviously  we  cannot  merely  set  w  =  ^, 
but  it  is  Dlausible  that  setting  ^  =  some  function  of  w 
might  do  the  trick,   VJhat  we  want  is  a  function  ^  =  ^(w) 
such  that  the  equation 

(3.8)  ^(w)  =  f(z,a) 
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has  a  power  series  solution 

(3.9)  w  =  a^Ca)  +  a^(a)z  +  a^^^"^^      "^  ••• 
where 

(3.10)  a^(a)  ^   0 

for  all  a  in  a  neighborhood,  |a|  <  A^,  of  a  =  0  ( including 
a  =  0).   Well,  if  w  =  w(z)  has  a  power  series  expansion 
such  as  (3.9 ),  then 

(3.11)  a  (a)  =  g 
■^      ^^  z=0 

From  (3.8)  we  have 

(3.12)  K.^  (w)dw  =  f,(z,a)d2  , 

whence  A*' 

f  (z,a)  r^^c^' 

(3.13)  a  (a)  =  ^ .  /"" 

^'(^^   z=0  (.^'qi^''' 

Since  f  vanishes  for  z  =  z,{a)  and  z.(a)  — >  0  as 
z  +        + 

a  — *.  0,  a,  (a)  would  seem  to  —>  0  as  a  — >0.   There  is, 
however,  at  least  one  way  in  which  this  undesired  occurrence 
can  be  nullified,   ^'e  can  pick  the  function  ^(w)  such  that 
^'  (w)  =  0  when  w  =  w_^,  where  vj  are  the  values  of  w  which 
correspond  to  z  =  z^  through  equation  (3.8),  i.e.  ^^(w^)=f (z^,a) 
This  means  that  the  function  ^' (w)  must  have  zeros  at  at 
least  two  prescribed  points;  hence  ^(w)  must  be  at  least  a 
cubic.   Since  two  values  of  its  derivative  are  prescribed, 
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we  pick  a  cubic  depending  on  two  parameters,  namely  the 
following: 

O.lU)    g(w,a,b)  =  b  +  3a^w  -  w^  . 

We  therefore  set 

(3.15)     g(w,a,b)  =  b  +  3a^w  -  w^  =  f(z,a)  . 

From  (3.15)  we  have 

dw  _   ^z(^'<^)      ^2^^'^^ 


(3.16) 


dz   g^(w,a,b)   3(a2.„2j 


and  therefore,  since 

(3.17)  B^  =  0   '  -  -  w  =  +  a  , 
we  should  pick  a  and  b  such  that 

(3.18)  g(+a,a,b)  =  b  +  2a^  =  f(z+,a) 

since  then  f  and  g  will  vanish  simultaneously, 
z       w  '^ 

Thus, 

(3.19)  /■  b  =  |[f(z^,a)  +  f(z^,a)] 
and        '^ 

(3.20)  V  a^  =  ^[f(z^,a)  -  f(z^,a)3   . 

Equation  (3.l5)  with  a  and  b  given  by  (3.19)  and  (3.20) 
defines  w  as  an  implicit  function,  w  =  w(z,a).   It  will  be 
proved  that  this  function  can  be  expanded  into  a  convergent 
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power  series  in  the  neighborhood  of  the  origin,  a  series 
which  will  be  valid  for  all  a  in  a  neighborhood  of  a  =  0, 
including  a  =  0,  and  which  is  such  that  the  coefficient 
of  z  is  different  from  zero  for  all  such  a.   This  last 
fact  means,  of  course,  that  the  inverse  function  z  =  z(w,a) 
is  also  a  power  series  in  w  good  for  all  a  in  a  neighborhood 
of  a  =  0. 

We  shall  give  a  precise  formulation  and  rigorous 
proof  of  this  key  theorem.  For  this  purpose,  it  is  more 
convenient  to  recast  the  formulation  of  the  theorem  a 
little  bit.   First  of  all  we  transform  the  origins  of  the 
J,  and  w  planes  to  one  of  the  zeros  of  f  and  g  ,  say  z_^ 
and  a,  i.e.  we  introduce  the  new  variables 

(3.21)  X  =  z  -  z^ 

; 

and 

(3.22)  V,.  y  =  w  -  a 
and  we  call 

(3.23)  f(z,a)  =  f (x  +  z^,a)  =  F(x,a) 
and 

(3.2i|)      ^  =  z_  -  z^   . 

Then,  by  (1.3)  -  (1.7),  we  have 

F'(0,a)  =  f (z^,a)  =  0 
(3.25) 

F>(C,a)  =  f'(z_,a)  =  0 


26, 

,  P"(0,a)  =  f"(z^,a)  ^0    if  a  ?^  0 

(3.26)  \      F"(4,a)  =  f"(z  ,a)  ^  0    if  a  ?^  0 

I 

l^  P"(0,0)  =  f"(0,0)  =  0 

.  F"'(0,a)  =  f'"  (z.,a)  f^   0  ,   (all  a) 

(3.27)  i 

'.  p'*'(C.a)  =  f"'  (z.,a)  f^   0   (all  a). 

Hence  we  can  write  (Jm*^^* 

(3.28)  f(z.a)  =  f(2^,aj+i(z-z^>^f"(z^,a)  +|(2.a^)3f*"(z^,a)  +  .. 
or 

(3.29)  F(x,a)  =  P(0,a)  +  |  x^  P"(0,a)  +  |x^p"'  (0,a)  +  ...  , 

where  F"(0,a)  — >  0  as  a  — ?.  0,  but  F   (0,a)  does  not. 
There  is  no  loss  in  generality  in  assuming  that  ^F"(0,a)  =  a, 
for  this  merely  amounts  to  replacing  a  by  the  new  parameter 
«P"(0,a),  a  replacement  v/hich  in  no  way  affects  any  of  the 
important  features  of  the  problem.   Then  calling 

(3.30)  P(0,a)  =  p, 

we  can  v.Tite  (3.29)  in  the  form 

(3.31)  F(x,a)  =  p  +  ax^  +  YX'^p(x,a)  ,    Y  ?^  0  , 

where  p(0,0)  =  1  and  p(x,a)  is  analytic  in  x  and  continuous 
in  a  in  a  neighborhood  of  x  =  0,  a  =  0, 

Correspondingly,  the  function  g(w,a,b),  defined  by 
(3.1i4-)  goes  over  under  the  change  of  variable  (3.22)  into 
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g(w,a,b)  =  (b  +  2a^)  -  3a  y^  -  y^  , 
or,  letting 

(3.32)  b*  =  b  +  2a^  , 
we  have 

(3.33)  g(w,a,b)  =  G(y,a,b'"')  =  b*  -  3a  y^  -  y^  . 

Here,  a  Is  the  same  as  the  a  in  g.   By  (3.20),  (3.23),  (3.2i^) 
and  (3.30)  this  means  that 

(3.31|)  a^  =^[^(2  +  ,^)  -  f(z.,a)l  , 

or 

(3.35)  a^  =  ^[F(0,a)  -  F(4,a)3  , 
or 

(3.36)  a^  =  |[p  -  F(?,al  , 

while  b'"'  is  obtained  by  using  (3.19),  (3.20),  (3.23)  and 
(3.32),  as 

(3.37)  b"  =  b  +  2a^  =  f(2^,a)  =  F(0,a)  =  p  . 

The  relations  (3.36)  and  (3.37)  giving  a  and  b'""  in 
terms  of  p  and  B,   could  also  have  been  obtained  directly 
by  requiring  that  G  and  F  vanish  at  corresponding  points, 

y     X 

for  from  (3.25)  we  have 

(3.38)  F'(x,a)  =  0 ',•■•■'">   X  =  0  or  X  =  4 
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(for  X  and  a  in  a  sufficiently  s'nall  neighborhood  of 
X  =  0,  a  =  0),  while  from  (3.33), 

(3.39)     G-(y,a,b''')  =  -3y(2a  +  y)  ; 

hence 

(3.i|.0)   G'(y,a,b""')  =  0  ^"-^  y  =  0  or  y  =  -2a. 

a  and  b"  are  then  to  be  picked  such  that 

(3.i;l)    G(0,a,b'')  =  F(0,a) 

and 

(3.I4.2)    G(-2a,a,b''-')  =  P(?,a)  , 

i.e.  such  that 

(3.i;3)    b'^  =  F(0,a)  =  p 

and 

(3.Ui|)    b"''"  -  i+a^  =  F(C,a)  , 

which  agree  v/ith  (3.36)  and  (3.37). 
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§U.   The  Key  Theorem 

Using  (3.23),  (3.31)  and  (3.33),  equation  (3.l5) 
takes  the  form 

(i;.l)   b"'  -  3ay^  -  y^  =  P  +  ax^  +  YX^p(x,a),  y  ?^  0  . 

But  b"  =  p,  (see  (3.37)),  and  therefore 

(1^.2)   -3ay^  -  y-^  =  ax^  +  YX^p(x,a)  ,   y  7^  0  , 

which  is  the  same  equation  as  would  be  obtained  ifp 
were  zero,  V/e  can  therefore  assume  p  =  0  to  begin  with 
and  in  addition  we  can  replace  F  by  -P  which  clearly  does 
not  alter  the  truth  or  falsity  of  the  theorem  we  wish  to 
prove,   •'•he  equation  to  be  treated  is  thus  reduced  to 

(I4..3)    3ay^  +  y-^  =  ax^  +  YX-^p(x,a)  ,   y  7^  0  , 
where,  since  p  =  0  and  F  has  been  replaced  by  -F, 

(k.k)       a  =  jF(4,a)  =  ^[a4^  +  Y4^p(C,a)]   (by  (3.36))  . 

Having  dispensed  with  the  preliminaries,  we  can  now 
formulate  the  main 
Theorem!   Let 

(l|.5)    P(x,a)  5  ax^  -t-  YX'^p(x,a)  ,  y  f^  0   , 

where  p(x,a)  is  analytic  in  x  and  continuous  in  a  in  a 
neighborhood  of  x  =  0,  a  =  0,  and  where 

(I;. 6)    p(0,0)  =  1  . 


liiji'v.; 
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Let 

ik'l)  G(y,a)  =  3ay^  +  y^  , 

and  let  4  be  the  root  of  F' (x,a)  =  0  (looked  at  as  an 
equation  In  x)  closest  to  but  distinct  for  a  ;^  0  from  the 
origin,  X  =  0.   (If  a  is  sufficiently  small,  such  a  root 
always  exists  and  is  unique  by  virtue  of  the  form  (i^..5) 
assumed  for  F(x,a), 

Then  if  and  only  if  a  is  chosen  such  that 

(1;.8)   a^  =  |F(^,a)  , 

the    equation 

(!;.9)        G(y,a)    =  F(x,a)    , 

i.e. 

(«4-.10)  3ay^  +  y^   =  ax^  +  Yx-^p(x,a)    ,      Y  7^  0    , 

has  a  unique  solution, 

(il..ll)    y  =  H(x,a) 

such  that  H(0,a)  =  0  for  all  a  in  a  neighborhood  of  a  =  0 
and  such  that  H(x,a)  is  a  power  series  in  x, 

(lj..l2)   y  =  H(x,a)  =  h^(a)x  +  h2(a)x  +  h^(a)x"^  +  ...  , 

good  for  X  and  a  In  a  neighborhood  of  x  =  0,  a  =  0,  a 
power  series  in  which 

(U.13)    h^(a)  ^  0 
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for  all  a  In  a  neighborhood  of  a  =  0.   Note  that  here,  as 
elsewhere  in  this  discussion,  we  employ  the  usual  definition 
of  nelghborhoodj  i.e.,  a  =  0  is  included  in  its  own  neigh- 
borhood.  If  we  wish  to  indicate  that  a  =  0  does  not  belong 
to  its  neighborhood,  we  shall  use  the  standard  terra  "deleted 
neighborhood"  of  a  =  0. 

The  inequality  (1^.13)  means,  of  course,  that  the 
inverse  function,  x  =  h"  (y,a->,  is  also  analytic  in  y  in 
a  neighborhood  of  y  =  0  and  has  a  power  series  expansion 
about  y  =  0  with  a  non-vanishing  coefficient  of  y  for  all 
sufficiently  small  a,   i'hus,  (I;. 11)  gives  a  one  to  one 
conformal  mapping  of  a  neighborhood  of  the  origin  in  the 
complex  x-plane  onto  a  neighborhood  of  the  origin  in  the 
complex  y-plane, 
P r oof  of  the  Theorem : 

We  have  already  seen  that  the  conditions  of  the 
Theorem  are  necessary.   The  real  problem  is  to  prove  their 
sufficiency.   To  do  this,  we  start  in  the  natural  way; 
we  use  ik'B)    in  (i|.,8).   This  gives 

(l+.ll;)   a^  =  ^[aE,^   +  Y4^p(4,a)3  . 

On  the  other  hand,  from  the  definition  of  4  (see 
p.  30)  we  have 

(k-lS)  2a  +   3Y4p(4,a)  +  Y4^P'(4,a)  =  0  . 
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If  the  value  of  4  in  terms  of  a  determined  by  (^.15) 
is  substituted  into  (i|..l'4),  then  a  is  determined  as  a 
function  of  a.   It  is  somewhat  more  convenient,  however, 
to  use  (I|..l5)  to  eliminate  a  from  (it..li4.)  and  look  upon  the 
resulting  equation  together  with  equation  (i;.l5)  as  para- 
metric equations  for  a  and  a  in  terms  of  the  parameter  ^. 
Note  that  a  — >  0  entails  E,   —oO. 

Using  the  abbreviations 

(U.16)    p(4,ct)  =  p   and  p»(4,a.)  =  p» 

(we  can  ignore  the  dependence  of  p  upon  a  because  p  does 
not  tend  to  zero  with  a)  we  obtain  from  (i;,l5) 

(U.17)   a  =  -  |y  4  P  -  |y  4^P'   , 
or 

(I+.18)   a  =  -  IyIOp  +  4p')  . 

Substituting  (i;.17)  into  (l|.,ll^.)  we  get 

(U.19)    a  =  -  |y^^^(p  +  Cp')  ^^  . 

Putting  (I;. 18)  and  (1;.19)  back  into  (I;. 10)  we  then 
obtain 

(1^.20)    -  |y  ^  4(P  +  4p'  )  ^3  y^  +  y^  =  -  |y4(3p  +  4p«  )x^ 

+  Ya  p(x,a)  * 
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In    (I+.20),   put 

(!;.21)        y  =  Y    -^  (P  +  Cp'  )     ^    x  +  v(x,a)    . 
Then 

-  |y(p   +  4p'  )4x^  -  3y  ^^  (P  +  4p'  )  ^Cxv 

-  |y  ^  (P   +  ^P' )  ^    ^v^   +  y(p    +  4p' )x^ 

+  3y  ^'^  (P  +  Sp'  )  ^    x^v  +  3y    -^  (P  +  ^P'  )         xv^ 

+  V^  =  -  |y(3p  +  4p' )4x^  +  YX^p(x,a)    , 
or 

([;.22)       3y^3(p  +  4pi  )  6  (x-4)xv  +  3y   -^  (p  +  4p' )     ^(x-|^)v^ 

+  v^  =  Y     \[p(x,a)    -  (p  +  4p')]x  +  4^P'^      x^   . 

But,    by   Taylor's   Theorem, 

(!^»23)        p(x,a)   =  p   +   (x  -  C)p'    +  (x  -  §)^  R(x,^)    . 

Therefore  the  coefficient  of  x   in  ([|..22)  can  be  written 
in  the  form 

Y  •  [(x  -  4)p'  +  (x  -  4)^  ^(x,C)  -  4p']x  +  4^p\ 


=  yCx  -  4)^[p'  +  xR(x,4)3   , 


,:\f    ;■• 
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so  that  (i;.22)  becomes 
{k'2k)        3y  ^^  (pHp'  )  ^^  (x-4)xv  +  3y  ^  (p+4p'  )  ^(x  -  |4)v^ 

+  v^  =  y(x-£J^  [p'  +  xR(x,4)lx^  . 
In  ik'^k)   put 
(i^..25)   V  =  x(x-i)u  . 
Then  Hi., 21^)   becomes 

(/4.26)   3y^(p+4P'  )  ^3u  +  3y  ^  (P+CP'  )  ^^  U  -   |fJu^ 

+  x(x-4)u^  =  Ytp'  +  xR(x,4)3  . 

Since  p  =  p(C,a)  — ->  1  as  a  — ?>  0,  the  coefficient  of 
the  lowest  order  term  in  u  in  (i4..26)  Is  bounded  away  from 
zero  as  a  --->0,  which,  by  the  remark  at  the  top  of  p.  32  is 
therefore  also  true  as  E,   — *►  0.   Since  equation  ([^..26)  is  of 
the  formn(u,x,a)  =  0  where i^  is  analytic  in  u  and  x  and 
continuous  in  a  and  where  i"^  (0, 0,a)  =  0,  this  means  that 
.  '  (0,0, a)  ^   0  for  all  sufficiently  small  a.   The  solution 
u  =  u(x,a)  to  such  an  equation  is  given  by  (see  [l] ) 

(         u  ^--^^/^^'^^  du 
(i+.27)   u  =  u(x,a)=^       7rr^^''^ 

|u|=e 

where  |u|=e  is  a  circle  sufficiently  small  to  ensure  that 
■,_(u,x,a)  7^  0  thereon  (since  i'i^(0,0,a)  j^  0,  p.(u,x,a)  ^  0 
in  a  neighborhood  of  u  =  0  and  such  a  circle  can  therefore 
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always  be  chosen. ) 

From  (i|,27)  it  follows  that  u  is  analytic  in  x  and 
continuous  in  a.   The  coefficients,  u  (a),  in  the  power  series 
expansion  for  u! 

CD 

([)..28)    u(x,a)  =  y-     n(a)x'^ 

m=0  " 

are  given  according  to  Cauchy's  Theorem  by 

„  i„\    _  ml    (     u(x,a)  , 
|x|=0 

and  are  therefore  continuous  in  a  because  u  is. 
From  ([|..25)  and  ('4.28)  it  follows  that 

(!;.29)   v(x,a)  =  ^  u  (a)x"'"'^(x-4)  , 

m=0 

and   therefore   by    (1^,21)    that 

(I1.30)        y  =   [y/3(p  +  4pi)  /3_  ^  UQ(a)]x 

+  x^[u_(a)    +     f^     u   (a)x'"'-^(x-4)]    . 
°  m  =  1     "^ 

Since  4  — >  0  and  p  — >  1  as  a  — >0,  it  is  clear  that 
the  coefficient  of  x  in  ('4.29)  is  different  from  zero  for 
all  a  sufficiently  small  and  surely  for  a  =  0.   i'hus ,  if 
we  take  H(x,a)  as  the  power  series  given  by  the  right  hand 
side  of  (i|.29)  we  see  that  this  H(x,a)  satisfies  all  the 
desired  conditions  and  the  proof  of  the  main  theorem  is 
therefore  comolete. 
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§5.   Derivation  and  Proof  of  the  Asymptotic  Expansion 

With  the  aid  of  the  foregoing  Theorem  we  are  now  in 
a  position  to  extend  the  sadc'le  point  method  to  the  case  of 
the  coalescence  of  two  branch  points.   We  start  out  exactly 
as  van  der  V'aerden  does  by  making  the  change  of  variable 
(1.2)  in  the  integral  (1.1).   As  in  §2,  this  gives  an 
integral  of  the  type  (2.2),   Next  assume  that  the  two 
branch  points  of  the  mapping  (1.2)  which  lie  closest  to 
the  curve  t  '  on  the  right  depend  upon  a  in  such  a  way  that 
they  coalesce  when  a  =  0;  call  these  two  points  ^^(a)  and 
their  pre  images  z_^(a).   Since  ^(a)  are  the  two  branch 
points  closest  to  '  on  the  right,  we  can  deform  \    into 
a  path  I  '  which  passes  through  ^^  and  by  bending  it  a 
little  if  necessary,  we  can  make  )  '  pass  through  <q=^_j^(0)=^__(0), 
the  point  of  coalescence,  too. 

^-plane 
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The  deformation  of  I  '  into  P   is  to  be  carried  out  in 
such  a  i^ay  that  no  branch  points  or  poles  of  the  integrand 
are  crossed  over  in  the  process.   As  remarked  above,  the 
properties  ^f  ^_^  and  ^_  ensure  that  this  can  always  be  done. 

Finally,  we  assume  that  the  z^(a)  satisfy  all  the 
conditions  listed  on  pp»  3-k»      Then  we  know  by  the  Theorem 
and  the  discussion  preceding  it  on  pp.  ?.l\.-29   that  if  we 
perform  the  additional  mapping 

(5.1)   g(w,a,b)  =  b  +  3a^w  -  w^  =  -?:  =  -f(z,a)  , 

then  z.  are  not  branch  points  of  the  mapping  from  the 
z-plane  to  the  w-plane.   All  other  branch  points  of  (1.2) 
remain  branch  points  under  (5.1),  however,  for  g^  =  0 
only  at  w  =  +  a,  and  these  points  correspond  to  z. .   At 
all  other  points  z  ^   z^  such  that  f  (z  , a)  =0,  g„  7^  0, 
and  such  points  therefore  remain  branch  points. 

Since  w  =  +  a  corresponds  to  z  =  z^  under  (5.1) 
and  since  by  hypothesis  z  —^z-^   and  a  — ^>  0  as  a  — >  0, 
of.  (l.ij.),  it  follows  that  w  =  0  corresponds  under  (5.1) 
to  z  =  Zq.   The  curve,  i"'    ,  in  the  ^-plane  will  therefore 
go  over  into  some  curve,  K,  in  the  w-plane,  which  will 
pass  through  the  points  w  =  ±  a  and  w  =  0  and  which  will 
once  again  go  to  infinity  in  both  directions.   In  fact, 
it  is  easy  to  determine  the  directions  in  which  K  goes 
to  infinity.  According  to  (5.1) 
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w  ->^  ■>        for  large  ^. 


Now  on  1  ,  ^  .-...,  I^je  '  in  one  direction  (namely  in  the 
upper  half-plane),  where  0  <  :,,.'<  i  and  ^       -^  |e  ^^^"  '  '  ^^ 
in  the  other  direction,  where   0  </.•  <  ^  , 
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Consequently,  if  w^|w|e^-^  in  one  direction  on 

in  the  other  direction  on  K,  then 


K,  and  w^-jwle 


whence 


and 


J-l  =  2  ;   XL'  =  21^ 


0  <  n  <  I 


f  <  n'  <  f 


-I  ir, 
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Therefore  K  is  a  curve  looking  something  like  the 
one  pictured  on  the  preceding  page.   The  integral  (1,1), 
meanwhile,  goes  over  into 

(5.2)  F(n,a)  =  e"^  j  e^^^^^^"^^^  4(w,a)dw  , 
where 

(5.3)  4(w,a)  =  ^[z(w,a)]^  . 

Again  we  need  to  impose  a  mild  restriction  such  as 

{2»l\.)»      In  the  present  case  an  even  weaker  assumption  than 

(2,i|)  suffices;  viz.  we  assume  that 

•J 
(5.ij.)    I4(w,a)|  <  m'  e^'^'     (m'  and  >v  constants) 

sufficiently  far  out  on  K.   This  being  the  case,  we  can 
deform  K  into  the  path  L_^  =  I^  -  I_^  shown  on  page  7, 
so  that 

f  2        3 

(5.5)   F(n,a)  =  e^^   I    e"^^^  ^"^  >  4(w,a)dw   . 

K 

Now  assume  that  a  (and  therewith  a)  is  sufficiently 
small  to  ensure  the  validity  of  the  Main  Theorem,   Then  z 
and  therewith  (1»  =  '^(^j^^'gn  i^  an  analytic  function  of  w, 
possessing  a  power  series  expansion  converging  in  some 
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circle   |w|  <  B  about  the  origin.   If  a  is  sufficiently 
small,  then  the  point  w  =  3a  and  therefore  surely  w  =  a 
lies  within  this  circle  (note  that  our  proof  of  the  ■'■'heorem 
showsthat  the  radius  of  the  circle  cannot  tend  to  zero  as 
a  — >  0;  this,  in  fact,  is  the  inoortant  point  of  the 
Theorem).   Hence  4  is  analytic  in  a  neighborhood  of  w  =  a. 
So  if  we  write 

(5.6)  (|>(w,a)  =  ^^{\]   ,a)  +  w  ^^(w  ,a)  , 

(which  can  always  be  done),  then  4-1  (■'.■»'i)  and  (j)p(''!»ci)  are 

2 
analytic  in  a  neighborhood  of  *;  =  a  and  consequently 

possess  power  series  expansions  of  the  form 

(5.7)  4i{'/,a)  =  ^  c^^^(a)(a2  -  ^,>)   ,   1  =  1,2, 

^  '      k=0   ^ 

whence 

(5.8)  4.(w^,a)  =  Y-     ci^^a)(a2-w2)   ,   i  =  1,2. 

As  explained  above,  the  series  (5.8)  converge  in  a  neigh- 
borhood  |w|  <  B  of  the  origin  including  the  points 
w  =  +a.   We  can  therefore  write 

K  k 

(5.9)  4.(v;^,a)  =  y"  ci^^a)(a2-w2)   +r,(w^,a)  ,  (1  =  1,2), 

where 


a^'s   i'^o'ijlis'-) 


k2. 
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(5.10)  rj^(w^,a)    =,A.(w,a)(a'^-w^) 


K+1 


where         lA(w,a)|   <  _/Vfor      |w|    <  B*. 

Next   take   a   circle      |w|    =  B  of  radius      |3a|    <   B  <   B* 

about   the  origin  and  denote   by  L^   the   portion  of  L 

If 
lying   in   this    circle   and  by  L_^  the   remainder   of  L_^.      Then 

we   can  write 


w-plane 


i+3. 


(5.11)  P(n,a)    =  e^^[F^(n,a)    +  F2(n,a)]    , 
where 

(5.12)  F^(n,a) 
and 


r 


J 


^nOa^w-w^)    ^(„^^)d^ 


(5.13)  F2(n,a)    =      \      e'^^^^^^"^    ^    (i(w,a)dw      . 


The  behavior   of  Fp  for   large   n  can  be   estimated   im- 

mediately.      On  the  part   I-  ©f  L^  running  along   Iq    (of. 
diagrams   on  p,    7   and  p,   \\2)    we   have 

w  =    C   ,    where  ^    is   real   and  0<Bg(l    <oo. 


On  the   part,    I_^,    of  L     along   I    , 
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w  =   e^^  L  ,   r' 


real   and  0  <  B  <   /l   <  oo  , 


Hence 


|F2(n,a)|   g 


/^' 
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„   e^^^^^^-"^)    4(w,a)dw 


I 


„   e^^^a^w-w^)   ^(^^^)^„ 
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r    ^n(3a2r-r^)    4(^^,a)dr 
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,n(3a2   e   ^    c  -  T^ )    ^(e"^r,a)dc 


But  since  i3a|  <  B, 
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(5.11+)      i3a^ti  <  ^B^iri  g  ^irP         for  irl  ^  B. 

2Tii 
Similarly,     |3a      e  -*     '^ I    <  tItI"    •      Furthermore,    by    (5.U) 


iS.k)  I4(w,a)|    <  M*    e^l^l 


3 


on  the  path,   (  C    is  real  and  positivej  so  we  can  omit  the 

3  K 
absolute  value  signs).   Therefore  if  n  >  ^-^  , 

,   f   (K-ln)r^ 
(5.15)    |F2(n,a)|  <  2M    \      e  ^  dT  . 

B 
If  B  ^  1,    then   it   follows   from   (5.15)    that 

^  (K-   2n)^3 

|F2(n,a)|    <   2M     j  3r^   e  ^  dT 

B 

Thus,    if  n  >  ^  +  1, 

3      -   ^^  B-^ 
|,F2(n,a)  |    <  2m'    e  ^^     e     "^ 

Calling 


and 


2m'    e  ^^^  =  M-^ 


2B^  .  ^ 


I  •>»- 


i: 


fe-'' 


we  have 

|P2(n,a)|  <  M*  e"^" 

■a  \jr 
for  n  >  ==^-p-   +  1,   In  other  words, 

(5.16)    P2(n,a)  =  0(e"'i'^).     (q  >  0) 

On  the  other  hand,  if  B  <  1,  then  we  write  {5«l5) 
in  the  form 


|P2(n,a)|<  2M 


e    ^  dC  ■»•[   e     ^  dl 

''b  1 


The  second  integral  here  is  0(e~^  )  by  the  argument 
just  given  above.   As  for  the  first,  we  have 

C^      (K-  |n)r^       (K-  |n)B3  f^ 

I   e     ^      dr^  e     ^         dr 

^B  "  -^B 


(>t-  ln)B^ 
=  (l-B)e     ^  =  OCe"^"^)  . 


Thus,  in  any  case,  (5,16)  holds. 

The  estimate  (5.16)  will  enable  us  to  show  that  the 
dominant  contribution  to  the  asjn^ptotic  expansion  of  F  is 
made  by  F, .   In  fact,  the  actual  expansion  is  obtained 
from  P,  by  inserting  the  series  expansions,  (5.9),  into 
(5.6)  and  the  result  into  (5.12).   This  gives 


I|6. 

(5.17)  F,(n,a)  =  Y:     4^^a)  f   (a2.w2)ke^<3a2w.w3)  ^^ 
•^        k=0   ^      ^  ' 

^+ 


L+ 


+  /,  r^(w^a)  e"(3^^^"^ 


^  dw 


L 

+ 


K  K 

=  n  4^^a)A(n,a)+  ^  4^^a)  Bj^(n,a) 


+  Rj^(n,a)  +  R2(n,a), 

In  order  to  transform  the  integrals  A,  and  B,  into 
Integrals  which  can  be  easily  evaluated,  write 

(5.18)   A;(n,a)  =  [^    (a^-w^)"  e"'^^^''-"^)  d„ 


L* 


■/ 


(a2-„2)''  ,n(3a^w-w3)  ^„  .  [  ^^^.„^^\n^i^^«.«h^„^ 


II 


Since 


vr 


".  I  J  :^  ; 


■Ttiy 


■, ; 


o 


■rfl 


,v, 


t.i 
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the  same  argument  as  used  on  F^Cnja)  (cf.  pp.  k3-kS) 
applies  to  the  second  integral  here  and  shows  that  it  is 
0(e"^'^).   Hence 

(5.19)  \in.a)   =  f   (a^-v,^)^  ^n(3a2w-w3)  ^^  ^  ^^^-qn^   ^ 

Similarly, 

(5.20)  B^(n,a)  =  j   w(a2.w2)^  ^nOa^w-w^)  ^^  ^  o{e-'*")  . 

So  if  we  call 

(5.21)  Aj^(n,a)  =  J   (a^.v,^)''  ^r^i^^^^'^^) 


+ 
and 


(5.22)   B.  (n,a)  =  f    ^,{a^.^^)^   ^nOa^w-w^)  ^^  ^ 


we  have 


L+ 


K 


(5.23)   Fj^(n,a)  =  f^  [cj^^^a)Aj^(n,a)  +  c^^^a)Bj^(n,a)] 


-qn. 


+  R^(n,a)  -f  R2(n,a)  +  0(e'^  )  , 
or,  with  the  further  abbreviation 


(5.21;)   Cj^(n,a)  =  c^^^a)Aj^(n,a)  +  cl^\{n,&) 


'  { 


•1  ./ 


■;( 


■■"-•-''      .r'-'* 


•■.  "^  i^  -'- ' 
>,'■■»*■'•• 


,  't  i 


;    f ,  .-^ . 


■  f     ■>  • . 


'p-  5  ^ 


'^f 


,>•',;       . 
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• "  »« 


(•-'•.' 


^vat: 
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K  _„ 

(5.25)  PT(n,a)  =  7"  Ci,(n,a)  +  R,(n,a)  +  Rp(n,a)  +  OCe""*"). 

i       k=0   **        ^  '^ 

We  wish  to  determine  the  order  of  the  terms  In  (5«25) 
for  large  n.   To  do  this  we  first  integrate  (5.21)  by  parts. 
If  k  ^  2,  we  obtain 

A^(n,a)  =  -  ^  /  (k-l)(.2w)(a2.w2)^"V(3aVw3)  ^^^ 

Thus, 

(5.26)  Aj^(n,a)  =  ^^^^  B^^^^^"'^^  • 

If  k  =  1,  then  A, (n,a)  =  0  because  then  the  integrand  is  a 
perfect  differential.   Hence,  if  we  agree  to  define  B_,  =  0, 

(5.26)  holds  for  k  ^  1. 

Similarly,  by  integrating  (5.22)  by  parts,  we  get 

(5.27)  B^(n,a)  =  -  ^  /  e'^(3aVw3)  ^^,(,2.^2^^-13  ,^ 

=  .^/e-(3a2w.w3)^^2.^2)k-\^, 


L+ 


+  2(|Z1)  /e'^(3aVw3)  ^ 


K 


^(a^-w^)    dw   (if  k  >  2) 


2 


=  -  h  *k-l<"'«'  -  ^TH^  Ak.i<n.»)  *  ^^TH^  A,.2(n,a), 


.  > ; 


U9. 
so  that 

(5.28)   Bj^(n,a)  =  ^^-^^^  Aj^.^(n,a)  ^  ^^^^^^   Aj^_2(n,a)  . 

This  holds  if  k  s  2.   If  k  =  1,  then  we  see  from 

(5.27)  that  B^  =  -  JL  Aq.   So  if  we  define  A^^^  =  0,  then 

(5.28)  holds  for  k  ^  1. 

Inserting  (5.28)  into  (5.26),  we  have 

(5.29)  A  (n,a)  =  ^^'^   [2{k-3)&\   .  (n,a )  -  (  2k-5)A.  ,(n,a)3, 
and,  on  the  other  hand,  if  we  put  (5.26)  into  (5.28)  we  get 

(5.30)  B.  (n,a)  =  2(k>l)a^  >  2(k-3).  3   ^^^^^ 

^  On)'^         ^"^ 

(2k-l)2(k-2)  n   /„  «\ 
_ B,  ,(n,a). 

In  other  words  A,  and  B,  both  satisfy  recursion  formulae 
of  the  type 

(5.31)  D^  =  ^  (XD^.^  *   ,xD^.3) 

where  here  and  in  what  follows,  \  and  |j,  denote  generic  co- 
efficients depending  on  k  but  not  on  n.   In  order  to  save 
writing  we  shall  call  all  such  coefficients  X  not  bothering 
to  distinguish  one  from  another  by  subscripts  or  the  like. 

Since  A,  and  B.  satisfy  a  relation  such  as  (5.31),  so 
too  must  C^   (see  (5.2I4.)), 

°k '  -^  '^S-u  *  ^S-3' 


>,Vj  a- 
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^  TlnT^  ^^^k-12  "■  ^3^^k-ll'^k-10»S-9^^' 
or,  in  general 

where  f^  is  linear  in  all  its  arguments. 


Now 


Hence 


A:,  =  0  ,   B,  =  ^  Aq  . 


^1  ■  3H^o  • 


^2=^%'      ^2-lk^O' 


Therefore,  ^2  "  ;:^  ^0  * 


Hence 


A-5  =  ^  B,  =  — ^5  A^  ; 


3   5n  2   3n  1    (^n)^ 


^3  ^  73^ '°  ■ 


Prom  (I+.32)  it  therefore  follows  that  if  k  =  l+m 


v-,r^'s' 


■>' I   -r 


•    ( 


;•)        i  :. 


■>  ■*■ 
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;    'm  I'i.j'-: 


'  ••  "  ,.•*  ;<      m  I't.i! 


'-. . 
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If  k  =  Um+1, 


Vn  =  -^  ^1  *  -~Si  4  ^  "(n-^^-'X^o  . 


or 


If  k  =  lim+2,   then 


(5.35)    0(^.2  =  -^5;^  C,  .  o(n-2"-i)Co  . 
If  k  =  l4in+3,   then 


^Uin+3 

(3n)^  ^^3  ^  (3„ 

j^+'2  ^0  ■" 

o(n' 

■2-2  )C, 

(5.36) 

^lpi+3  ' 

=  ,3„,2-2  S  ^  o< 

-2m-2  vp 
n             )Cq 

• 

Collecting 

;  (5.33)  -  (5.36), 

vje  have 

< 

1 

_  k 

0(n     2)Cq 

if 

k  = 

km 

(5.37) 

k+1 

l^x  = 

0(n        2    )Cq 

if 

k  = 

km+1 

:^^^ 


J  I 
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C.  =  0(n   ^)C.  If  k  =  Um+2 


'k   "^"    '"0 
(5.37) 


or 


k+1 
^  Cj^  =  0(n  ~^)Cq  if  k  =  Uin+3  , 


^k  ~  ^^"  ^^^0  ^^  ^  ^^  even. 


k+1 


(5.37') 

It  remains  to  estimate  the  remainder, 
(5.38)   R,(n,a)  =  J   e"*^''^"-"^'  r,  (w^a)dw  . 

To  do  this  we  use  (5.10).   Since  4  Is  analytic  and  of 
exponential  order  at  oo  ,  so  too  is./L-,   Just  as  in  previous 
cases  it  therefore  follows  that 


(5.39)   Ri(n,a)  =  j   e'^^^^  """  ^J^.^{^^,a){&^•^,^}        dw 

L 


+  0(e""^)       (q  >  0). 

For  simplicity,  assume  that  K+1  is  even,   (If  K+1  is 
odd,  only  minor  changes  in  the  following  argument  are 
necessary, ) 


V   ■  - 


•^O 


'-!'-         > 


,)■  . 
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,   Ci. 
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.■;    i^:  Jj^ 
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Integrating  (5.39)  by  parts  gives 

Now  it  is  clear  that  we  can  continue  to  integrate  by 

2  2 
parts  until  we  reach  an  integral  in  which  (a  -w  )  no  longer 

appears,  since  each  integration  by  parts  decreases  by  one 

2  2 
the  highest  power  of  (a  -w  )  appearing.   It  is  also  clear 

that  this  will  happen  first  when  we  apply  all  the  ^'s 

occurring  in  the  repeated  integration  by  parts  to  the  power 

2  2 
of  (a  -w  )  left.   If  we  apply  some  of  the  differentiations 

to  the  -A.-J ,  then  we  shall  have  to  carry  out  more  integrations 

2   2 
by  parts  in  order  to  eliminate  the  (a  -w  )  from  the  integrand, 

These  additional  integrations  by  parts  will  increase  the 

pov;er  of  —  multiplying  the  resulting  integrals.   Since  K+1 

is  assumed  even,  it  is  therefore  evident  from  the  treatment 

of  the  C.  integrals  that 

^  K+1    /      2   3 

(5.1+1)   Ri(n,a)  =  0(n  "^)  J      e""^^^   """^  ^TV],  (v;,a)dw  . 

Once  again  we  use  the  fact  that  the  part  of  the 
integral  taken  over  V^   is  0(e~"'^).   This  means  that 

-  ^   f       2   3 
R^(n,a)  =  0(n   ^    )  j      e"^^^  ^"^  U^(w,a)dw  . 


.  t-ii.'-^    ^(6    i  \'  .   ' 


'  1       vih 


n. 


L.    *■ 


WO" 


.)v.'     .■ 


.  ■'   l-ii. 


»*.-:        ..I'i. 


:j.'' 


H  rJ    »■•  Ci 
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.   ».      ,_  , 


iAd.i    aim 


n 
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But 


o'i  L^,     lAi(w*a)  I   <  _A_«      Hence 


K+1 


lRl(n,a)|   ^  0(n       ^   )  j         \ 


^nOaVw^)^^,      ^ 


and   therefore,    sxarely 


_  K+1       /^  2       3 

(5.1^2)  |R;L(n,a)|   ^  0(n     "^)^      |e"^^^  ^"^    ^ 

L. 


dw      , 


Similarly,    for 


RAT\,a)   =  /       e 


nOa^w-w-')..  _    ,..2 


w  r2(w   ,a )dw   , 


we  have 


(5.1+3)         |R2(n,a)|   g  0{n'      ^    )  j      j^g'^Oa  w-w    )^^| 


K+1        ,^ 
2 


,2,.  „3 


The  estimates  (5.^2)  and  (5.1+3)  differ  slightly  from 
those  obtained  for  the  individual  termsof  the  series. 
Since,  however. 


|AQ(n,a) I  = 


e"(3aVw3)^^ 


|Ao(w,a)|g  i   le^^^aVw^)^^,, 


,    !•><. 
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likewise 


|Bo(n,a)|  <  I   Iwe'^^^^^^-'^^^dwl   , 


K 


therefore,  if  we  put 


;j;=  J      le^naVw^)^^,  ^  ^  j^^nOaVw^) 


C.. 


dw|  , 


+ 


then  obviously 


K+1 


(S.kk)  |Ri(n,a)|  g  0(n   ^  )Cq 

and  (cf.  (5.21;)) 

(5.i;5)     |CQ(n,a)|  ^  m"Cq   . 

Relations  (5.37'),  iS.kk)   and  (5.U5)  establish  the 
asymptotic  character  of  the  series  (S»2S)*      This,  in  turn, 
together  with  (5.11)  and  (5.16),  establishes  the  validity 
of  (5.23)  as  an  asymptotic  representation  for  F(n,a)  and 
completes  the  theoretical  part  of  the  discussion. 
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§6.   Calculation  of  the  Terma  in  the  Lxpanslon 

We  turn  next  to  the  practical  side  of  tlie  problem 
and  the  first  question  which  naturally  arises  is  how  to 
compute  the  terms  in  our  asymptotic  expansion  (1.1!;). 
What  have  to  be  computed  are  the  c5   (a)  and  the  A.  's  and 
Bj^'s.   According  to  (5.7) 

(6.1)  J^>(a)=i:i)l!cfl^)(a2,a) 

'^  kl    ^ 

where  ^.      (a  ,a)  denotes  the  k   derivative  of  4^  v;ith 

2 
respect  to  its  first  argument  evaluated  at  a  . 

For  the  purpose  of  the  applications  it  Is  somewhat 

more  convenient  to  express  the  c)    In  terms  of  4  instead 

of  4,  and  c|p»   Since  4-i  and  w^p  are  just  the  even  and  odd 

parts  of  4>  this  is  not  hard: 

(6.2)  4i(w2,a)  =  |f4(w,a)  +  4(-w,a)] 

(6.3)  <L(w^,a)  =  ^[cj.(w,a)  -  4(-w,a)]  . 
Hence 

(6.i|)    c^^^a)  =  4i(a^,a)  =  |[4(a,a)  +  4(-a.a)l 

(6^5)    cj,^^(a)  =  <^^(&^,&)   =  -^[(^(a.a)  -  4(-a,a)j   . 

The  cl^     with  k  >  0  can  be  computed  from  (6.2)  and  (6.3)  by 


■i;: 
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2 
differentiating  with  respect  to  w  and  then  setting  w  =  a. 

It  should  be  noted,  however,  that  in  the  computation 

of  the  coefficients  c)   (a)  by  means  of  formulas  (6.2)-(6.5) 

indeterminate  expressions  such  as  rr  will  always  occur. 

These  then  have  to  be  evaluated  by  using  L'Hopital's  rule 

once  if  a  ;^  0  and  twice  if  a  =  0.   That  this  is  necessary 

should  be  no  surprise,  for  it  will  be  recalled  that  at 

the  very  outset  when  we  looked  for  a  method  to  resolve  the 

difficulty  of  the  coaloscence  of  two  branch  points  one  of 

the  key  ideas  was  to  make  the  numerator  as  v;ell  as  the 

dz 

denominator  of  the  expression  for  -t-   vanish  and  thereby 

render  the  expression  indeterminate  in  the  hope  that  it 
would  then  approach  a  finite  non-zero  limit  when  a  — ?>  0. 
That  this  hope  can  be  realized  has  been  proved  above  and 
we  therefore  know  that  the  coefficients  can  be  computed  by 
L'Hopital's  rule  and  will  all  be  finite  with  the  first  one 
of  them  different  from  zero  (some  after  the  first  might  be 
zero,  however,  which  is  no  contradiction  and  no  harm). 

As  an  illustration  of  these  ideas,  we  shall  compute 
ci   (a),  the  first  and  most  important  coefficient.   According 
to  (1,12) 


(6.6)    4(+a,a)  =  l(z+)f| 


w=+a 
z=z . 


We  assume  for  simplicity  that  J  is  regular  and  ?^  0  at  z^. 
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For  the  application  we  have  in  mind  ^  is  in  fact  a  constant, 
We  therefore  only  have  to  compute 


dz 
3w 


w=+a 
z=z . 


lim  ^ 


Z-e^Z 


11.  l^lz-I) 


Z-*»Z 


(by  (1.8)) 


Thus, 


(6.7)    lim 
w-s»+a 

z-^  z . 


lim 


6w 


dz 


w-^ta  f,,(z,a)-^ 


zz 


dw 


z— ^z 


dz 
dw 


+6  a 


(by  L'Hopital's  rule). 


for  a  ^  0, 


^zzz^V""^ 


Therefore, 

(6.8)   4(ia,a)  =  f(z^) 

So,  by  (6.14.), 


+6a 


if  a  =  0. 


(^.^)  4'^-hi^^.^/T^^7^^U^^j/r;pj^ 


;    n  i     ..  ■:.   *>■    v!.n .       < . 


>  i^'' 
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It  remains  to  express  the  A,  »s  and  B,  's  in  terms  of 
some  well  known  functions.   Since,  as  remarked  In  the 
introduction,  the  recursion  formulas  (1.17)  and  (1.18) 
(or  (5.26)  and  (5.28))  permit  all  the  A^»s  and  B^'s  with 


k  a  1  to  be  expressed  in  terms  of  A^  and  Bq,  it  is  only 
necessary  to  express  A_ 
This  is  done  by  putting 


necessary  to  express  A-  and  Bq  in  terms  of  known  functions, 


(6.10) 


-^3 
w  =  (3n)  and 


a  =  (3nr^^^/2 


in  (1.15)  (or  (5.21)).   Then 

-^3  f     zt-  k3 
(6.11)   AQ(n,a)  =  (3n)    J     e         ^   dt 


or 


(6.12)    AQ(n,a)  =  ^  Ai^(2)   , 

(3n)^3 

where  Ai_j_(z)  is  the  Airy  integral  defined  by  (1.20).   Thus, 

(3n) 


(6.13)    A^(n,a)  =   ^"^^ ^  -  Ai^[(3n)  ^^   a^]  ,  which  is  (1.21) 
0  1/3    + 


As  for  Bq, 


(6.11+)    B^ 

L 


o(n,a)  =  J      we"(3^^^-^^>  dw   (by  (1.22)) 


+ 


5L't 


(pi 
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-  -jpj  -—^   J       e         dw 


oa   ^ 


so  that  by  (6.13) 

(6.15)  B.(n,a)  =  _2Si   Al'[(3n)  "^^  a^]  , 

(3n?5 

which  ia  (1.22) 

To  obtain  (1.21;),  merely  put 

lit 

(6.16)  t  =  e  "^r' 


in  (1.20).   Then 


(6.17)   Ai(z)  =  %^i   j   e^  3-   dr, 


in   Z'  lit 

T  f   e"T  zr  -.  H^ 


whence  by  (1,23), 

(6.18)  Ai(z)  =  ^  2k       ^1^^  "^  z)» 
so  that 

-  ^  -i^ 

(6.19)  AQ(n,a)  =  ^^  ^  .    h-L(e  "^  z)  , 

k(3n)  ^3 
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or 


lit  lit  2 

(6.20)  AQ(n,a)   =     ^   ^'^ ^    h^[e"  "^(3n)^  a^]    , 

which   is    (1.21;). 

From   (6.1U)    and    (6.20),    we    then  get    (1.25): 

'^  .171  2 

(6.21)  BQ(n,a)  =  ^  ^     .     -  h^[e     "^(3n)^  a^]    . 

k(3n) '^ 

As  mentioned  in  the  introduction,  formulas  (1.2i|)  and 
(1.25)  are  the  more  useful  ones  for  actual  numerical  compu- 
tations.  Also  important  is  the  observation  that  although 

3 

only  Q.      is  specified  by  (1.9),  there  is  no  latitude  in  the 

3  - 

choice  of  a,  because-  /a  must  be  chosen  so  as  to  be  con- 
sistent v/ith  the  choice  of  cube  roots  made  on  pages  37-39. 
In  ge^sral,  it  is  best  to  e::anine  each  individula  case  on 
its  own  merits,  taking  into  account  the  particular  function 
f(z,a)  involved.   This  point  vill  be  further  illustrated 
in  §8. 
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S7»   Parameters  Not  Independent 

The  formula  (l.l!;)  and  the  special  case  of  it  (1.28) 
are  good  for  all  sufficiently  small  values  of  a  (and  there- 
with a)  and  sufficiently  large  values  of  n  and,  as  remarked 
many  times  before,  retain  their  validity  even  as  a  -*  0, 
They  were  derived  under  the  assumption  that  a  and  n  were 
independent  parameters.   If,  however,  a  depends  upon  n  in 
such  a  way  that  a  —e>  0  as  n  — «►  oo  ,  then  it  is  clear  that 
(l.lij.)  and  (1.28)  are  also  applicable  in  this  case  since  it 
is  but  a  special  case  of  the  more  general  one  v;hich  has 
been  treated  above.   For  example,  if 

(7.1)  a  =  a(n) 
where 

(7.2)  lim   a(n)  =  0, 
n  — 5»CD 

then  a,  being  a  function  of  a,  is  also  a  function  of  n: 

(7.3)  a  =  a(n) 
such  that 


(7.i|.)     lim   a(n)  =  0  , 
n — ^>CD 


and  (1,28)  becomes 


'.nvl   9flT 
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(7.5)   nnMnn^^.-±-^)iu^^/;=^r^ 


2k(3n)^3 


zz 


_    iTt  2 

e     "^   (3n)^  a2(n) 


If  a(n)   happens   to  be   a   function  of  n  such  that 


(7.6) 


a(n)   =  0(n     ^)    , 


then,  to  lowest  order,  the  argument  of  h-,  is  independent 
of  n  and  the  asymptotic  representation  simplifies  somewhat. 
By  making  an  assumption  which  is  more  or  less  equivalent 
to  (7.6),  F,  W,  J,  Olver  [8]  has  obtained  asymptotic 
expansions  for  Bessel  functions  of  nearly  equal  order  and 
argument  which  are  sl-nilar  to  those  obtained  in  the  next 
section  by  applying  the  method  developed  above. 
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§8.  Application  to  Hankel  Functlona 
The  Hankel  functions  [2], 


(8.1) 


Hil'(z)  = 


1  I  ^-l(z  sine  -  n9)  ^,  ^ 


where  C  is  the  path  shown  below i 


A 


9-plane 


-It 


<    C     < 


.  v: 


<^t} 


{en 


U»8) 


65. 


have  been  the  subject  of  numerous  investigations,  a  number 
of  which  have  had  as  their  object  the  determination  of  the 
asymptotic  behavior  of  H    (z).   It  is  found  that  the 
principal  cases  of  interest  are  those  In  which  both  the 
order,  n,  and  the  argument,  z,  become  large  and,  in  fact, 
at  the  same  rate.   More  precisely,  the  important  cases  are 
those  in  which  [11] 


(8.2)   H^l>(z)  =  H^l)(nr)  =  -  i  /  e'^^^^  ^ ^^  "  ^^  d© 


where  r  is  a  fixed  positive  number  independent  of  n.   This 
is  of  the  form  (1.1)  with  z  called  9,  a  called  r, 

(8.3)    f(^,r)  =  i(r  sine  -  ©) 
and 

{Q.k)        f(e)  =  -  ^  . 

From  (8.3)  we  have 

(8.5)  f'(e,r)  =  i(r  COS0  -  1)  . 
Consequently,  the  branch  points  of  the  mapping 

(8.6)  K,   =  -f(e,r) 

are  given  by  the  roots  of  tire  equation 

(8.7)  cos©  =  i  . 
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To  determine  their  order,  we  have  to  look  at 

(8.8)  f"(©,r)  =  -Ir  sinG  . 

Prom  (8.7)  and  (8.8)  it  Is  clear  that  all  branch 
points  are  simple  if  r  ;^  1  but  at  least  double  if  r  =  1. 
In  fact,  since 

(8.9)  f"*(©,r)  =  -ir  cos©  , 

they  are  precisely  double  if  r  =  1, 

We  therefore  have  a  situation  in  which  all  branch 
points  change  their  order  as  the  parameter,  in  this  case 
r,  tends  to  a  certain  value  (one).   VJe  have  already  seen, 
however,  that  if  we  perform  the  mapoing  (8.6),  then  the 
dominant  terms  in  the  asymptotic  expansion  of  H   (nr) 
will  be  contributed  by  those  branch  points  in  the  9-plane 
whose  images  in  the  ^-plane  lie  closest  on  the  right  to 
the  curve  P  in  the  ^-plane  which  is  the  image  of  C  under 
(8.6).   It  is  well  known  from  the  theory  of  Bessel  functions 
that  it  is  the  branch  point  ©  =  Q_  closest  to  the  origin 
from  which  the  main  part  of  the  asymptotic  expansion  arises. 
According  to  the  discussion  of  sections  1-6  it  will  there- 
fore suffice  to  remove  the  point  ©^  as  a  branch  point. 
Naturally  the  branch  point  -©,  with  which  it  coalesces  must 
also  be  removed.   We  therefore  take  for  the  points  @_^ 

(8.10)  ©^  =  +6^ 


e,^., 
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where  e^^  Is  the  root  of  (8.7)  having  s-nallest  absolute 
value. 

Since 

(8.11)  r  coseQ  =  1  , 

it  follows  from  (8.3)  that 

(8.12)  f(e^,r)  =  +i(r  sinS^  -  9^)   =   +i(tan  ©^  -  9^) 

=  ±i  (^7^1   -   tan-^  /T:!') 


/l  -  ^  dp   . 
■^1 

Hence,  according  to  (1.9)  and  (1.10) 

(8.13)   a^  =  -  |(tan©Q  -  9^)  =  -  ^(^7^  -   tan"^  ^7^') 

and 

(S.lij.)    b  =  0, 

while  according  to  (8.8)  and  (8.12) 


(8.15)   f"(teQ,r)  =  :;:i  /r^-l   . 


To  determine  v/hlch  cube  root  to  take  in  (8.13),  note 
that  for  a  =  0,  (1.8)  becomes  (cf.(8.3)) 

-W 7- ■»-..., 


or 


w  =  (|)"^3  Q   + 


-^iiJiv-: 
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For  a  ^  0,  y ^^ 

3a^w  -  w^  =  l(r-l)9  +  ...  ,  .^^ 

-  r  I'' 

w  =  i(i^)e  +  ...  .  -  4*X 

3a^ 

We  therefore  have  to  choose  the  cube  root  in  a  such  that 

i  ^  — >   (r)  ^  (where  we  take  1  ^=  e^  ), 
3a        ° 

P  V^ 
i.e.,  since  r  is  real,  such  that   a  i  -^  is  pure  imaginary 

and  has  imaginary  part  of  the  same  sign  as  r-1.   This  means 
Ini 

if  r  >  1,  and 

271  i 

a  =  '^ 


^-t7-  (/T^*  -  tanh"Vl-r^r3 

if  r  <  1,   Therefore,  formula  (1.28)  tells  us  that  if  r  >  1, 
then 

(l)/-^^  -i  /-6(/r^-l'  -   tan'VT^)    ^ 

V-5    /  1/ 

k(3n)/3v/  2/3^^2:^ 

where 


/fl  i/L\    u(l)^„^^            "i            /-D(/r   -1   -   tan    /r   -1)    -^         v.    /    n 
(0,16)   H        {T\r}^ =:— ~    / — ^ i: .    h,  (z), 

"  V3   v/  ■^'  -^ 


.    (3n)'/3     (./?:r  ■  tan-VTIT)  ^3 


2 

If  r  <  1,    then  , 

(8.17)  H'l^nr),^— i^^/ii^SZ^^anl^VSTLl  .   h^,,,, 

k(3n)7       a'/^/x:^' 

where 

_   (3n)'^3(/i,r^'  -   tanh'Vl-r^')  ^3 


.8rt 


:,''K'/ 
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(8.18) 


For  r  =  1,   we   have  to  use    (l.?7).      This   gives 

V3    • 


n 


k(3n) 


As  mentioned  previously,  the  ordinary  saddle  point  method 
gives  rise  to  three  different  formulae,  each  valid  for  a  dif- 
ferent range  of  r.   These  9re  the  following  (cf.,  'Vatson  [11] 
(pp.  S43-4)): 


(8.19)        H^^^nr)--® 


in(\/r^-l'  -   tan"Vr^,-l')  -  iiri 


J^nnjr^^l 


(r   >   1); 


(8.20) 


H^^W)/~- 


nC/l-r^   -  tanlT'^Jl^) 


v'2-irnJl-r^ 


^gn(tanh"-^Jl-r^   -  v/l-r^) 


/gTTn\ll-r^ 


(r   <   1); 


(8.21) 


nh 


"n      ("^^_2/3^„l/3 
2       Tin  ' 


1      _  ^31/3 


Formula  (8.2l)  is,  in  fact,  identical  with  (8.18),  as  is 
readily  seen  by  using  the  definition  (I.23)  of  hj^(z)  in  (8.18) 

The  following  table  and  graphs  offer  a  comparison 
among  the  actual  values  of  K^      (nr),  those  coriputed  by 
using  the  uniform  expansion  ( (8.16) -(8.18 ) )  and  those 
computed  by  using  the  ordinary  expansion  ( (8.19 )-(8,21) ) : 
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H^^^^nr)  =  J^(nr)  +iY^(nr) 

r 

__  1 

n 

Actual  Values 

Values  Using 
Uniform  Expansion 

Values  Using 
Ordinary  Expansions 

2 

3 

.1147+. 3282  1 

.1144+. 3310  1 

.1040+. 3342  i 

4 

-.1053+. 2829  1 

-.1058+. 2818  1 

-.1133+. 2812  i 

5 

-.2340+. 1354  1 

-.2337+. 1345  i 

-.2376+. 1306  1 

10 

.1864-. 0438  1 

.1856-. 0435  1 

.1871-. 0419  1 

1.5 

3 

.4247-. 0090  1 

.4240-. 0074  1 

.4352+. 0199  i 

4 

.3576+. 0983  1 

.3570+. 0991  i 

.3578+. 1196  i 

5 

.2834+. 1754  1 

.2826+. 1756  i 

.2786+. 1904  i 

10 

-.0900+. 2199  1 

-.0882+. 2162  1 

-.0959+. 2185  1 

1.1 

3 

.3587-. 4291  1 

.3595-. 4273  i 

.5206-, 4382  1 

4 

.3364-. 3698  1 

.3369-. 3688  i 

.4615-. 3665  1 

5 

.3209-. 3260  1 

.3211-. 3253  i 

.4220-. 3159  i 

10 

.2804-. 1983  i 

.2803-. 1979  1 

.3271-. 1786  i 

1.0 

3 

.3090-. 5385  1 

.3101-. 5371  1 

.3101-. 5371  1 

4 

.2811-. 4889  1 

.2817-. 4880  1 

.2817-. 4880  1 

5 

.2611-. 4536  i 

.2615-. 4530  i 

.2615-. 4530  i 

10 

.2074-. 3598  1 

.2076-. 3596  i 

.2076-. 3596  i 

.95 

3 

.2818-. 5969  1 

.2754-. 6093  1 

.3923-. 8371  1 

4 

.2507-. 5522  i 

.2561-. 5433  i 

.3361-. 7328  i 

5 

.2280-. 5214  1 

.2426-. 4955  1 

.2974-. 6623  i 

10 

.1661-. 4439  1 

.1789-. 4177  i 

.1992-. 4945  i 

.  f 
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H(^)(nr) 

r 

n 

Actual  Values 

Values  Using 
Uniform  Expansion 

Values  Using 
Ordinary  "LXpansions 

.9 

3 

.2540-, 6600  1 

.2643-.  6i|.l5i 

.3068-. 7402  i 

4 

.2197-. 6215  i 

.2118-. 6386  1 

.2575-. 6614  i 

5 

.1947-. 5963  i 

.2061-. 5733  i 

.2233-. 6104  1 

10 

.1246-. 5454  i 

.1133-. 5798  i 

.1350-. 5046  i 

.75 

3 

.1711-. 9101  1 

.1591-. 6219  i 

.1728-. 7719  i 

4 

.1320-. 9167  1 

.1370-. 8970  1 

,1309-. 7643  1 

5 

.1046-. 9445  1 

.0641-. 9634  i 

.1024-. 7816  i 

10 

.0421-1.2769  1 

.0348-1.0191  i 

.0370-1.0798  i 

— p 
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Actual    value 
Uniform    expansion 


Ordinary   expansion 


V^^) 
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1.0 


— I — 
1.4 
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— I — 
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-I— 
8 


1.0 


— I — 
1.2 


— I — 
1.4 


— I — 
1.6 


— r- 
1.8 


2.0    ^ 


|Y3(3r) 


.2  - 
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f  Y4(4r) 


1.6  1.8      ''       2.0 


78. 


tY.^dOr) 
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2.0 
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